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Abstract 

The telegraph process models a random motion with finite velocity and 
it is usually proposed as an alternative to diffusion models. The process 
describes the position of a particle moving on the real line, alternatively with 
constant velocity +v or —v. The changes of direction are governed by an 
homogeneous Poisson process with rate A > 0. In this paper, we consider 
a change point estimation problem for the rate of the underlying Poisson 
process by means of least squares method. The consistency and the rate of 
convergence for the change point estimator are obtained and its asymptotic 
distribution is derived. Applications to real data are also presented. 
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1 Introduction 



The telegraph process describes a random motion with finite velocity and it is usu- 
ally proposed as an alternative to classical diffusion models (see Goldstein, 1951 
and Kac, 1974). The process defines the position of a particle initially located at 
the origin of the real line and moving alternatively with constant velocity +v or 
—V. The changes of direction are governed by an homogeneous Poisson process 
with rate A > 0. The telegraph process or telegrapher's process is defined as 



where V{0) is the initial velocity taking values ±v with equal probability and 
independently of the Poisson process {N(t),t > 0}. Many authors analyzed 
probabilistic properties of the process over the years (see for example Orsingher, 
1990, 1995; Foong and Kanno, 1994; Stadje and Zacks, 2004; Zacks 2004). 
Di Crescenzo and Pellerey (2002) proposed the geometric telegraph process as 
a model to describe the dynamics of the price of risky assets where X{t) replaces 
the standard Brownian motion of the original Black-Scholes-Merton model. Con- 
versely to the geometric Brownian motion, given that X(t) is of bounded varia- 
tion, so is the geometric telegraph process. This seems a realistic way to model 
paths of assets in the financial markets. Mazza and RuUiere (2004) linked the 
process (11.11) and the ruin processes in the context of risk theory. Di Masi et al. 
(1994) proposed to model the volatility of financial markets in terms of the tele- 
graph process. Ratanov (2004, 2005) proposed to model financial markets using 
a telegraph process with two intensities \± and two velocities v±. The telegraph 
process has also been used in ecology to model population dynamics (see Holmes 
et ah, 1994) and the displacement of wild animals on the soil. In particular, this 
model is chosen because it preserves the property of animals to move at finite ve- 
locity and for a certain period along one direction (see e.g. Holmes, 1993, for an 
account). 

For the telegraph process {X{t), < t < T} observed at equidistant discrete 
times Q = < . . . < tn, with tj = zA„, z = 0, . . . , n, riA^ = T and A„ ^ as 
77, oo, De Gregorio and lacus (2006) proposed pseudo-maximum likelihood and 
implicit moment based estimators for the rate A of the telegraph process. Under 
the additional condition nA„ — > oo as n — > oo, lacus and Yoshida (2007) studied 
the asymptotic properties of explicit moment type estimators and further propose 
a consistent, asymptotically gaussian and asymptotically efficient estimator based 
on the increments of the process. 

In this paper we suppose that for a telegraph process occurs a switch of the rate 
from Al to A2 at some time instant G [0, T] and the interest is in the estimation 
of the change point 6^0 and both Ai and A2. 




t > 0, 



(1.1) 
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The change point estimation theory has been employed widely by means of 
the likelihood function (see Csorgo and Horvath, 1997). Unfortunately, the like- 
lihood function for the telegraph process is not known and the pseudo likelihood 
proposed in De Gregorio and lacus (2006) is not easy to treat in this framework. 
We will then proceed using the alternative method based on least squares proposed 
in Bai (1994, 1997) and used in different contexts by many authors including Hsu 
(1977, 1979) for the i.i.d. case and Chen et al. (2005) for the mixing case. Our 
model is peculiar in itself for the properties of X{t) described in the above, be- 
cause A is a parameter related to the mean and the variance of the Poisson process 
and because the mesh A„ plays a role in the definition of the rate of convergence 
of our estimators. 

The paper is organized as follows. Section [2] describes the model, the obser- 
vation scheme and the change point estimator. The consistency of change point 
estimator is discussed in Section [3] and distributional results are presented in Sec- 
tion m Finally, Section [5] contains an application of our method to real data: we 
will reanalyze the classical change point data sets of Dow-Jones weekly closing 
(see Hsu, 1977, 1979) and IBM stock prices (see Box and Jenkins, 1970). 

2 Observation scheme and least squares estimator 

We consider a telegraph process {X(t),Q < t < T < oo} defined as in (11.11) 
and assume to observe its trajectory only in ri + 1 equidistant discrete times 
< ti < . . . < t„, with ti = iAn, i = 1, . . . ,n and nAn = T. We assume that 
a rate of changes of direction shift occurs during the interval [0, T] at an unknown 
time 6^0 = rT, r G (0, 1). Therefore the changes of direction are governed by 
an inhomogeneous Poisson process with parameter X{t) = Xil{t<eo} + ^2^{t>eo} 
where the positive values Ai, A2 and the change point r (or ^^o) are unknown and 
to be estimated given the observations X(ti), X(t2), . . . , X{tn). In order to sim- 
plify the formulas we use the following notation: X{ti) = X(iA„) = Xj. The 
asymptotic framework is the following: A„ — > and nA„ = T^ooasr7,^oo. 

The telegraph process is not Markovian and, as mentioned in the Introduction, 
it is not possible to derive the explicit likelihood function of the observations X/s, 
therefore we can not apply the statistical methods based on the likelihood function. 
To work out our estimation problem, we shall follow the approach developed in 
Bai (1994), which involves least squares-type estimators. The same point of view 
has been applied by Chen et al, 2005, in a context of financial time series. For our 
model, the time increment A„ plays an active role in the study of the asymptotics 
of our estimators so the proofs, although in some cases along the lines of Bai 
(1994) require some technical, but crucial, adjustments. 
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In order to obtain our estimator we introduce some basic notations. Let 



Y. - '^{\m\<v^n} . _ -, 

where rji = X-i — Xi^i is the increment between two consecutive observations. We 
indicate the mean value of Yi with 7^ = ^-e^'" = m = 1, 2. We 

observe that the random variables Yi are independent and identically distributed 
because depend by the increments |77j|. lacus and Yoshida (2007) proved that the 
estimators 

1 -^1 
In^-^Yi and A„ = -— log(l - 7„A„) 

are consistent, gaussian and asymptotically efficient estimators of 7 and A respec- 
tively. We will use these properties in the following without necessarily mention- 
ing them. 

We assume that the change occurs exactly at time U — tk^ — koAn — 9o, 
therefore ko — [nr], r G (0, 1), where [ ■ ] represents the integer- valued function. 
The change point estimator is obtained as follows 



^0 = argmin min <^ J^iYi - 71)^ + (^i - 72)' 

k \ 71.72 ^ — ' ^ — ' 

\ I i=l i=k+l 

( k n ^ 

= argmin \ ^(y, - Y^f + J] (y. - %^^f I , (2.1) 



1=1 i=k+\ 



where 



^^Y.{Y^-l.f = \Y.Y, = Y,, 
1=1 i=i 

n 1 " 

min V {Yi - 72)2 = V y^ = %-k- 

72 ^-^ n — k ^-^ 

i=fc+l j=fc+l 

We indicate the sum of the squares of residuals in the following manner 

k n 



u! = Y.^y^ - Ykf + E - (2.2) 



1=1 i=k+l 

then 



ko — arg min C/^ , (2.3) 

k 
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and 7i = F^^, 72 = Yn-ko respectively the least squares estimators of 71 and 
72- This gives the two estimators 



Ai = -— log(l -7iA, 



A2 = ---log(l-72A„). (2.4) 



By setting = ^ XlLi = ZlLi simple algebra leads to 



where 



and 



i=l 



2 l-'n-fe — Ik) 



\/k{n — k) 



SnDk 



^ k Sk 

J-'k — TT- 

n S„. 



Therefore, formula (I2.3h implies that 

ko = argmax V"^^ = argmax |-Dfe| = argmax \/k{n — k)\Vk 



(2.5) 



(2.6) 



(2.7) 



Our first result concerns the asymptotic distribution of ko under the null hy- 
pothesis that Ai = A2. This permits us to test if a shift has taken place during the 
interval [0,T]. 

Theorem 2.1. Under Hq, i.e. X2 = \i = Xq, we have that for — ^ 0,nA„ 00 
as n ^ 00 the following result hods 



nAr 



An 



I A 



\B'it)\, 



(2.8) 



where {_B°(t), < t < 1} zi' a Brownian bridge. 

Proof Let = - 7o, 7o = = ''"aT" = Aq + o(l). Then, = and 
al = Var{ii) 



— '^^-"-^^ ^ — i*— 5 — ° "-'^ — '5-^. We introduce the following function 



A2 



with 5„ = X]r=i ^i- '^^'■^ 



+ {nt - [nt]) — < t < 1, 



Vt 



a„.\/n 



i=l 



a. 



[nt] 



(2.9) 
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and 



Var 



[nt] . / 



Since 



= Var 
= t 

l{|„,|<.A„}-(l-e^O^")| 



Vi 



V 



'^{\Vi\<vA,^} 



Vint] An 



(2.10) 



/a: 



< 



the Lindeberg condition is true 



["*] p r -] ^ f 2 \ 

E -^ i-^{VA;rig,i>£V^iA;ro-n}^i / „ 

Then from (|Z9] ). (12.101) and (|2.11|) we can conclude that 

^ ^iV(0,t). 



(2.11) 



a„\/n 



(2.12) 



Now, by applying Donsker's theorem (invariance principle) we are able to 
write that 

X„(t) ^ B{t) 

{x4t)-tx4i)}SB\t), 

with B(t) and -B°(t) representing respectively a standard Brownian motion and a 
Brownian bridge. Let A; = [nt], we can write 



Xn{t)-tXn{l) 



1 



anVn 
1 

GnVn 



S 

Ok On 

n 



nt — [nt] 

H —C,[nt]+1 



crr,.\/n 



k 



1=1 



i=l 



nt — [nt] 
GnVn 



We observe that 



- 7o) - - Y.^Y^ - 7o) = - A 



and consequently 

\fnK'n\Dk\ 



V" Y- 



AT 



Xn(t) — tXn{l) LI^[„t]+l 



o-„.\/n 



(2.13) 



It is easy to see (by Chebyshev inequality) that 

nt — [nt] 



sup 

t 



0. 



By the law of large number ^ Y^^=i ~^ ^o, while 
Therefore from (12.131) follows that 



nAr 



An 



\Dk\^\B\t)\. 



(l_e^O'^»i)e-^0'^ii 



□ 



Corollary 2.1. The same convergence result of the Theorem 12. 1 [ follows when we 
consider 



VnKr, 



\Dk\, 



where Aq is any consistent estimator for Xq. 

Remark 2.1. From Theorem \2.1\ we derive immediately that for 6 G (0, 1/2) 



nAr 



A 



sup I Di, 



nAr 



5n<k<{l~5)n 
d 



sup |5°(t)|, 

S<t<(l-5) 



A, 



sup iVfcl A sup (t(l -t))-^/2|5°(t)|. 

Sn<k<{l-S)n 5<t<{l-S) 



(2.14) 



(2.15) 



The last asymptotic results are useful to test if doesn't exist a change point. In 
particular it is possible to obtain the asymptotic critical values for the distribution 
(|2.15l) by means of the same arguments used in Csorgo and Horvdth ( 1997), pag. 
25. 



3 The consistency properties of the estimator 

We shall study the consistency and the rate of convergence of the change point 
estimator (12. 7h . It is convenient to note that the rate of convergence is particularly 
important not only to describe how fast the estimator converges to the true value, 
but also to get the limiting distribution. The next Theorem represents our first 
result on the consistency. 

Theorem 3.1. The estimator f = ^ satisfies 

\f-T\ = {nAn)-^'\^2 - 7i)"'0p(v/l^) (3.1) 
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Proof. By the same arguments of Bai (1994), Section 3 and by using the formulas 
(10)-(14) therein, we have that 



|t - r| < a(72 - liY^ sup \Vk - EVk 

k 



(3.2) 



where Cr is a constant depending only on r. Let = "'"^'''''^"^"Z \ given 

V 

that 



k 1 



y/nAn V n — /c 



E2. 

i=A:+l 



we obtain that 



{Zn-k + Zk] 



(3.3) 



where Z^ = ^ Yli=i ^^d Z^-k = y^=^ EILfc+i ^y applying Hajek- 
Renyi inequality for martingales we have that 





Ell 


> a j 


, max 


l<fc<n 


Ck 





1 " 



E(Z,Y 



< 



AA„ + o(A„) >A 1 



fc=i 



A + 0(1) 1 

" fc=i 



(3.4) 



Choosing = Vk and observing that Efc=i k ^ < C log for some C > (see 
e.g. Bai, 1994), we have that 



max 



1 ^ 



n 



(3.5) 
□ 

Remark 3.1. By means of the law of iterated logarithm we obtain immediately 
the following rate of convergence which improve the previous result. We have that 



* 2=1 

Then from the relationships (13.31) and (13.51) we obtain the result (13.11) . 



T — T 



(nA„)-i/2^72 - 7i)"'Op(v/loglogn). 



(3.6) 
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Remark 3.2. Theorem \3.1\ implies that, under the additional hypothesis A„ = 
0(?7.^), e G (—1, 0) we have also consistency, i.e. (nA„)^(f — r) ^ in proba- 
bility for any j3 G (0, 1/2). 

We are able to improve the rate of convergence of f. 

Theorem 3.2. We have the following result 

1 



T-r = Or, 



nA„(72 - 7i)^ 



(3.7) 



Proof. We use the same framework of the proof of the Proposition 3 in Bai (1994), 
Section 4, therefore we omit the details. 

We choose a 5 > such that r G {5,1 — 5). Since k/n is consistent for r, 
for every e > 0, Pr{k/n ^ (5,1 — 5)} < e when n is large. In order to prove 
(|3.7I) it is sufficient to show that Pr{\T — r| > M{n^nlri)^^} small when n 
and M are large, where 7„ = 72 — 71. We are interested to study the behavior 
of Vk for n5 < k < n{l - 5), < 5 < 1. We define for any M > the set 
D^^M = {k : n5 < k < n{l - 5),\k - ko\ > MA"'^-/-^}. Then we have that 



Pr{\r-T\> M{nA^Ynr}<e + Pr{ sup l^^l > I^^J}, 



keD„ 



for every e > 0. Thus we study the behavior of Pr{snpk(zD^ \ Vk\ > \Vko\}. It is 
possible to prove that 



Pr { sup l^fcl > \Vk,\ \ < Pr'> sup Vk -Vk,>0 



+ Pr { sup Vk + Vk,<Q 



P + Q 



(3.8) 



Furthermore 



Q < 2Pr< sup 



fc<n(l-<5) n — k 



j=/c+l 



+2Pr \ sup — 

fc>n(5 



i=l 



By observing that Yl^m ^ ^ 



(3.9) 



0(m ^), the Hajek-Renyi inequality yields 



P { max 

k>m 



i=l 



> a y < 



1 (1 -e^^")e-^^'^ 



/a: 



(3.10) 
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where r.v.'s Zi are defined in the proof of Theorem 13. 1[ The inequality (13.101) im- 
plies that (13.91 ) tends to zero as n tends to infinity. Let h{k) = a/ {{k/n){l — k/n)), k 
1, 2, n, for the first term in the right-hand of (13.81) we have that 



P < Pr < sup 
-\- Pr I sup 



\H{k)\ > 



(3.11) 



Pi + P2, 



where 



i=l 



H 



^ n ^ n 



n — kn 



(3.12) 



(3.13) 



=fc+l i=ko+l 

We prove that Pi tends to zero when n and M are large. Thus we consider only 
k < ko or more precisely those values of k such that n5 < k < nr — MA^^7^^. 
For k > nS,we have 



kg - k 
n6ko 



^-^ n no ^-^ no ^-^ 

i=\ 2=1 i=k+l 

(3.14) 

where B >0 satisfies |6(A;o) - b{k)\ < B\ko - k\/n. By means of (1231), (ITTOl) 
and (|3.14l) . we obtain 



Pi < Pr { — 

nr 



[tit] 



i=l 



> 



+Pr { sup — 

i<k<n n 



+Pr 



1 



i=l 



> 



6B 



sup 



-2 nr — k 



< 



36D 36B'^D 
+ 



[nr] 

•t=fc+i 



> 



5'JnCr 



+ 



3QD 



where D 



(l-e 



. < A + o(l). When n and M are large the last three 

terms are negligible. Analogously we derive the proof of P2. 

□ 
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4 Asymptotic distributions 



We want to study in this Section the asymptotic distribution of f under our limiting 
framework for small variations of the rate of change of the direction. The case 
An = A2 — Ai equal to a constant is less interesting because when A„ is large the 
estimate of ko is quite precise. 

We note that A^ = A2 — Ai — implies 7n = 72 — 7i — ^ 0. By adding the 
condition 

A„^0, V nh ^ 

the consistency of f follows immediately either from Theorem 13.11 or Theorem 
13.21 In order to obtain the main result of this Section, it is useful to observe that 



ko = argmax\4^ = argmax?T,A„(\4^ — Vj^^) (4.2) 
and to define a two-sided Brownian motion W{v) in the following manner 

1^(.)4^^J7^' (4.3) 

where Wi, W2 are two independent Brownian motions. Now we present the fol- 
lowing convergence in distribution result. 

Theorem 4.1. Under assumption (|4.1I) . for nA„ — 00, A„ as n ^ 00, we 

have that 



A ° . t 2 _ 

where W{y) is a two-sided Brownian motion and A is any consistent estimator for 
Ai or A2. 

Proof. The proof follows the same steps in Bai (1994), Theorem 1, hence we 
only sketch the parts of the proof that differ. We consider only f < because of 
symmetry. Let Kn{v) = {k:k=[ko + v/\7^'^-i-\ -M < < 0, M > 0} and 

K{v) = nA^{V^ - Vl;) (4.5) 

with k G Kn{v). We note that 

nA„(lf-\4^J = 2nA„EVfc„(Vfe- Vfc„) 

+2n^r.{Vk,-EVu,){Vu-Vk,) 

+nAr,{Vk - \4o)' (4.6) 
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The last two terms in (14.61) are negligible on Kn{v). Since ^/nK^{ykQ — EVko) 
bounded by (13.21) . we have to show that A„ | Vk — Vk^ \ is bounded. In particular, 
we can write 

\/^n\Vk - Vk,\ < ^f^a\G{k) + H{k)\ + ^f^r.\EVk - EVk^ 

where G{k) and H{k) are defined respectively in (13.121) and (13.131) . The upper 
bound (13.141) is Op(l), because the first term is such that 



\pnKn 



ko — k 
n6kn 



ko 
i=l 



< 



M 



dm Anil 



ko 



J2 \^niY, - 



7i, 



i=l 



(4.7) 



similarly for the second term and for the third term we apply the invariance prin- 
ciple (12.121) . Now we explicit the limiting distribution for 

2nAnEVk,{Vk - Vk,) = 2^r{l - r)nA„7n(V[,„+,A-A-] - ^^o)- (4.8) 

For simplicity we shall assume that ko + t>A,~^7~^ and f A,~^7^^ are integers. We 
observe that 

nAniniVk - Vk,) = nAnin{G{k) + H{k)) - nAnin{EVk, - EVk), (4.9) 

where G{k), H{k) are defined in the expressions (I3.12|) . (|3.13l) . We can rewrite 
G{k) as foUws 



k-kQ ^ b{ko) - b{k) ^ 
-JkT - 7i) + 1 - 7i) 



G{k)= b{ko) 



1=1 



i=l 



ko 

+ b{ko)- $^(F,-7i). 



(4.10) 



=fc+i 



By the same arguments used to prove (14.71) we can show that the first two terms 
in (14.101) multiplied by nA„7„ are negligible on Kn{M). Furthermore 6(/co) = 



12 



V(l — r) and n/k ^ l/r for k G Kn{M), then we get that 

fco 

k 



i=fc+l 

I i=k+l } 

7„v/a;; Y ^n{Yi+k - 7i) > + Op{l] 



i=l 



^ 'XiWii-v) (4.11) 



where in the last step we have used the invariance principle (12.121) . Analogously 
we show that 



nA„7„i7(fco + t^A;S-2) ^ ^l^^ ^^w,{-v) . (4.12) 

i — r 



Since 



nAnin{EVu, - EVu) ^ , ' , (4.13) 



we obtain that 



In the same way, for f > 0, we can prove that 



K{v) A 2^-^\^W2{v) - Y| • (4.15) 

By applying the continuous mapping theorem and Theorem 13. 2[ 

nAn-il{f-T) d 1 



A Ai 



argmax A„(t>). (4.16) 



Since aW{v) = W{a'^v), a G M, a change in variable transforms argmax^, A„(f ) 
into Ai argmaxj; ) — , which concludes the proof. □ 

Using the consistency result, we are able to obtain the asymptotic distributions 
for the estimators Ai, A2, defined in (|2.4I) . 
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Theorem 4.2. Under the assumption (14.11) we have that 



where 



M ^'V(O.S), (4.17) 



Proof. We start noticing that 



^nA„(7i(A;)-7i(fco)) (4.19) 



1=1 1=1 



I h h 1 




i=ko 



Since fco = [™]: = fco + Op(A„-'^7^^), and r;,A„,7^ ^ oo, we have that 
(|4.19l) is (v^7n)~^0p(l), which converges to zero in probability. Then Ai(fc) = 
log(l -7i(A;)A„) and Ai(fco) = log(l - 7i(A;o)A„) have the same lim- 
iting distribution. Obviously the same result holds for A2. By Theorem 4.1 in 
lacus and Yoshida (2007), the convergence result (14.171) follows. 

□ 



5 Application to real data 

In this section we consider an application of our model to two well known real 
data sets. The first data set is about the Dow-Jones industrial average and the 
second one is the IBM stock prices. In both cases, the data mesh A„ is not close 
to zero, hence the asymptotics of our set up does not hold. Nevertheless, our 
findings seems to confirm the results of previous analyses. 
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Dow-Jones closings 
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Figure 1: Weekly closings of the Dow- Jones industrial average July 1971 - Aug 
1974 (up) and correspoind returns (down). 

5.1 Dow- Jones data 

This data set contains the weekly closings of the Dow- Jones industrial average in 
the period July 1971 - Aug 1974. These data have been proposed by Hsu (1977, 
1979) and used by many other authors to test change point estimators. There are 
162 data and the main evidence found by several authors is that a change in the 
variance occurred at point 89th which corresponds to the third week of March 
1973. Instead of working on the values we transform the data into returns as usual 
X{ti) = {W{ti) - W{ti-i))/W{ti^i), i = 1,.. . ,n with W the series of Dow- 
Jones closings and X the returns. We assume that X follows a telegraph process. 
In this application, the data are not sampled at high frequency, i.e. A„ is not close 
to zero, hence we test our estimator of the change point even if the asymptotics 
is not realized. Further, and for the same reason, we cannot assume as known the 
velocity v of the process hence we first estimate v by the average of the rescaled 
increments, i.e. 

t=i 

This is a consistent estimator of v, hence we construct the estimator of 7 as follows 

^^i^y; where y. ^ 

i=i 
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Dow-Jones returns 
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Figure 2: Change point estimates on the returns of the weekly closings of the 
Dow- Jones industrial average July 1971 - Aug 1974. Major change point estimate 
ko = 89 which corresponds to the 3rd week of March 1973; second change point 
estimate ki = 27, August 1971. 

With these quantities, we construct the statistics in (12.61) and maximize it. The 
maximum is reached at ko = 89 which confirms the evidence in Hsu (1974, 1979). 
Once we obtained the estimation of the change point, we re-estimate the velocity 
in both part of the series (before and after point 89th) and the two lambda's. We 
obtained respectively Vi = 0.61, Ai = 48.53 and V2 = 1.24 and A2 = 34.61 which 
confirms the intuition from the graphical inspection of the returns (i.e. in the first 
period there is a high number of switches but with low velocity which correspond 
to low variance of the returns; conversely for the second period). Looking better at 
the first part of the series, we observe that variance is not stable, so we re-run the 
procedure and obtained a new change point ki = 27 around august 1971. Figure 
[2] contains the two change point estimates plotted against the Dow- Jones returns. 

5.2 IBM stock prices 

This data set contains 369 closing stock prices of the IBM as, set. They have been 
analyzed in Box and Jenkins (1970) and further by Wichern et al. (1976) in order 
to discover change points. Box and Jenkins (1970) fitted an ARIMA(0,1,1) on the 
first order difference and discover heteroschedasticity; Wichern et al. (1976) fitted 
an AR(1) model on the first differences of the logarithms. We consider instead the 
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Figure 3: Return of the IBM stock closings (see e.g. Box and Jenkins, 1970). 
The major change point occurs at index ko = 235, the other two at ki = 18 and 
fca = 309. 

returns as in previous example and apply the same sequential procedure. Data 
are reported in Figure [3] along with a couple of change points discovered by our 
estimates. The first change point was found at point /cq = 235 which confirms the 
findings of Wichern et al. (1976). We further discovered another change point at 
time index ki = 18 on the time series on the left to ko and a second change point 
on the right-hand series at time k2 = 309. 
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